Precision Estimation of Crop Coefficient for Maize Cultivation Using High-Resolution Satellite Imagery to Enhance Evapotranspiration Assessment in Agriculture

The estimation of crop evapotranspiration (ETc) is crucial for irrigation water management, especially in arid regions. This can be particularly relevant in the Po Valley (Italy), where arable lands suffer from drought damages on an annual basis, causing drastic crop yield losses. This study presents a novel approach for vegetation-based estimation of crop evapotranspiration (ETc) for maize. Three years of high-resolution multispectral satellite (Sentinel-2)-based Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Normalized Difference Red Edge Index (NDRE), and Leaf Area Index (LAI) time series data were used to derive crop coefficients of maize in nine plots at the Acqua Campus experimental farm of Irrigation Consortium for the Emilia Romagna Canal (CER), Italy. Since certain vegetation indices (VIs) (such as NDVI) have an exponential nature compared to the other indices, both linear and power regression models were evaluated to estimate the crop coefficient (Kc). In the context of linear regression, the correlations between Food and Agriculture Organization (FAO)-based Kc and NDWI, NDRE, NDVI, and LAI-based Kc were 0.833, 0.870, 0.886, and 0.771, respectively. Strong correlation values in the case of power regression (NDWI: 0.876, NDRE: 0.872, NDVI: 0.888, LAI: 0.746) indicated an alternative approach to provide crop coefficients for the vegetation period. The VI-based ETc values were calculated using reference evapotranspiration (ET0) and VI-based Kc. The weather station data of CER were used to calculate ET0 based on Penman-Monteith estimation. Out of the Vis, NDWI and NDVI-based ETc performed the best both in the cases of linear (NDWI RMSE: 0.43 ± 0.12; NDVI RMSE: 0.43 ± 0.095) and power (NDWI RMSE: 0.44 ± 0.116; NDVI RMSE: 0.44 ± 0.103) approaches. The findings affirm the efficacy of the developed methodology in accurately assessing the evapotranspiration rate. Consequently, it offers a more refined temporal estimation of water requirements for maize cultivation in the region.


Introduction
The ongoing trends in climate change are exerting profound implications on various facets of human society, particularly impacting agriculture, and the global food supply [1].Improving agricultural water use efficiency and crop water productivity are the pillars of attaining higher agricultural production.Consequently, the optimal consumption of limited water resources is becoming of primary importance for the public sector, especially in South Europe [2].This is particularly true in the Mediterranean region, where multiple stakeholders (agriculture, household, industry) compete for increasingly scarce water resources [3].This is especially true in agricultural regions (e.g., northern Italy, Emilia Romagna), with water-intensive (maize, tomato) agricultural production.Changes in temperature, rainfall patterns, atmospheric carbon dioxide (CO 2 ), and the frequency and intensity of extreme weather events all have a substantial influence on agricultural yields, leading to rising food insecurity [4].Proper water management in agriculture is a key factor in agricultural production.
Due to the high spatiotemporal and temporal variability of meteorological parameters and vegetation as well as other environmental factors, such as soil conditions or the depth of the water table, planning optimal distribution of irrigation water is a very complex task [5].Such highly heterogeneous environments can be surveyed as detailed as possible resulting in a large variety of data.Meanwhile, data-intensive modern artificial intelligence, pattern recognition, and neural networks may contribute to more efficient planning and assessment of irrigation water distribution [6].On the other hand, the proper estimation of evapotranspiration (ETc) is a key parameter in the estimation of crop water demand for irrigation planning, design, and efficient irrigation water management of irrigation schemes [7,8].One of the most commonly used methods for ETc estimation is the crop coefficient approach, which can be calculated based on reference evapotranspiration (ET 0 ) and crop coefficient (K c ) [9].The K c changes throughout the growing season and reflects the proportion of the ground cover of the vegetation [10].Several approaches can be used to measure and/or estimate crop evapotranspiration (ETc), such as eddy covariance, soil water balance approach, Bowen ratio energy balance system, pan evaporation, and weighting lysimeters [11].However, these methods as well as ETc-based soil water balance models are expensive, hard to automate, and only provide point-based measurements.Hence, the application of these models encounters challenges when extended to expansive regions due to the spatial and temporal heterogeneity of input data.This is particularly notable as local conditions, including soil hydraulic parameters, interactions with groundwater, and variations in plant growth, pose intricate considerations that are not easily accounted for [12].
To tackle such challenges, remote sensing (RS) offers low-cost labor-saving, and timesaving approaches for assessing the spatial and temporal patterns of vegetation at different scales [13,14].Satellite-based global estimates of evapotranspiration are available like MOD16 (MODIS Global Evapotranspiration) with 1 km spatial resolution including both ET 0 , ETc estimates [15] or ECOsystem Spaceborne Thermal Radiometer Experiment on the Space Station (ECOSTRESS) with 70 m spatial resolution [16].The combination of different satellite data sources enables state-of-the-art and very sophisticated assessment methods.Traditionally, optical and thermal cameras (i.e., Surface Energy Balance Algorithm for Land (SEBAL) [17], Two Source Energy Balance (TSEB) [18], Surface Energy Balance System (SEBS) [19], Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) [20], and the Simplified Surface Energy Balance (SSEB) [21] and others) that are coupled with synthetic aperture radars can provide a variety of near real-time snapshots of the environment, which can contribute to modern agricultural methods [22][23][24].Earth observation (EO) approaches offer accurate data for precision agriculture applications and end-users (i.e., farmers, landowners, and decision-makers).High-resolution multispectral images may be used to schedule irrigation depending on crop demands in near realtime [25,26].Despite substantial research on ET crop estimates for water management using EO data, the poor geographical and temporal resolution of the sensors has hindered their usefulness and technical transfer [27].However, precision irrigation requires highresolution estimates of crop evapotranspiration and irrigation water requirements at the field level.There are several other hydrological models which use ETc as an input parameter as well [28].Using biomass-based daily estimation of K c over the vegetation period fosters a more reliable estimation of evapotranspiration contributing to the calculation of a more accurate water balance [29,30].In the crop coefficient approach, K c is typically taken from the Food and Agricultural Organization, FAO [31], which does not capture the spatial and temporal variability in estimating crop evapotranspiration.However, proper detection of crop growth affects the estimates of evapotranspiration fluxes and irrigation needs.Nowadays, high-resolution satellite datasets (10 m of spatial resolution) such as Sentinel-2 are easy and free to access online, therefore the spatial and temporal variability in the K c can be easily derived from VIs considering actual biomass, crop growth conditions, and stages [32].The Sentinel-2 A/B was used in this study, with five days of returning time.Free satellite imagery, high-resolution, innovative spectrum capabilities, a 290 km swath width, and frequent return durations are boosting operational and commercial uses of EO data for precision agriculture applications and research initiatives [33,34].The recent Sentinel-2 mission from the European Space Agency (ESA), as part of the Copernicus program [35], has significantly enhanced regular monitoring of agricultural parameters, like Normalized Difference Vegetation Index (NDVI), Normalized Difference Red Edge Index (NDRE), Normalized Difference Water Index (NDWI) and Leaf Area Index (LAI) [36][37][38].NDVI is a widely used VI and could provide important information for crop modeling research [39].NDVI time series are applied to estimate crop growth under different climatic conditions [40].The NDVI data has been used to monitor crop conditions in irrigation management [41] and estimate crop yield under different agroecological conditions [42][43][44].The spatial and temporal variability in K c can also be easily derived from NDVI because actual crop growth conditions and stages influence the temporal changes in NDVI values [32].The NDRE uses the reflectance from red edge bands to estimate plant health [45].The NDRE responds to plant health, biomass, and photosynthetic activity.It is, therefore, appropriate for plant health and productivity monitoring [46,47].The NDWI is a commonly used VI to determine the difference between water and vegetation cover, using the differences in the visible and near-infrared wavelength ranges [48].It is used to delineate, measure, and track the size and the spatial extent of change in water bodies [49] as well as to assess the effect of soil moisture or drought stress on plants [50].The monitoring of LAI distribution and variation is essential for tracking and detecting vegetation vitality and growth [51].LAI value may also be determined from RS data by empirical relationships between canopy reflectance and VI regardless of time or costs [52].
This research evaluates the use of RS-based indices in ETc calculation for the region of Emilia Romagna, Italy.There is an IT Irrigation Advisory Service for Farm Water Management services (IRRIFRAME) which utilizes a water balance model that aims at supporting crop irrigation management on a field scale and results in seasonal water savings of about 120 million m 3 in Italy [53].The IRRIFRAME provides evapotranspiration data under standard conditions [54], calculating reference evapotranspiration data (ET 0 ) by the Hargreaves formula with a 6.25 km 2 geographical grid.This study aimed to develop an advanced RS-based crop coefficient (K c ) and crop evapotranspiration (ETc) estimation method to support the preciseness of the calculation of the potential maize water consumption with 10 m spatial resolution.Sentinel-2-derived remotely sensed indices (i.e., NDVI, NDWI, NDRE, and LAI) were integrated with climatic data to define (i) the lengths of the phenological phases, (ii) biomass-based crop coefficient and (iii) to calculate and evaluate more accurate RS-based crop evapotranspiration.Due to the 5-day revisit time of Sentinel-2 satellites, the results contribute to the timely calculation of crop evapotranspiration.Results could also contribute to further improvement in spatial resolution (10 m) of the water balance calculations at maize fields.

Study Site
Field trials were carried out by the Irrigation Consortium for the Emilia Romagna Canal (CER) at the experimental farm Acqua Campus, covering 12.5 hectares of area that can be divided into approximately 25 fields.The fields are located in the plain of the Po River, in the province of Bologna, near the village of Mezzolara di Budrio (44   32 ′ E).Various crops are cultivated on the farm (both perennial and annual).Crop rotation in the experimental fields is characterized by annual crops, such as winter wheat, soya, maize, onion, and processing tomato.
The farm's soil is typical of the plain of the Po River and has a high content of silt, clay, and fine sand.The soil in the area can be described as clay/clayey loam.It belongs to the Italian soil group "Suoli SECCHIA franco argillosi", which can be classified as "Oxyaquic Haplustepts" (fine loamy, mixed, superactive, and mesic according to Soil Taxonomy).These soils are calcareous and moderately alkaline; they have a texture of clay loam in the superficial layers and loam deeper in the profile.Soil layers are affected by multiple floods of the nearby Idice River.The soil within the farm is heterogeneous, and its hydraulic characteristics vary from field to field and layer to layer, with a slight gradient from east to west.The mean values for the soil parameters are specified in Table 1.An extended, shallow groundwater table usually is present at a depth ranging from −0.6 to −1.8 m.During winter and at the beginning of the growing season, the capillary rise could be significant in terms of replenishing the evapotranspiration of crops.For the experimental farm, the water is supplied from the CER canal system, and its quality is regularly checked during the irrigation season [55].CER's irrigation water characteristics (long-term mean) are freely available online [56].
The farm has a fully operational weather station equipped with a rain gauge, anemometer, phreatimeter, and pan evaporimeter.
The climate of the site can be defined as sub-humid, with a mean annual temperature of 13.7 • C and a mean annual rainfall of 771 mm.The year 2018 of the field trials was characterized by severe drought.According to the estimates based on CER analyses, 2018 was the second harshest year in the study area, with a cumulative precipitation mostly below 200 mm in all areas which represents a deviation of 50% from the reference climate (period 1961-2018).These conditions contributed significantly to maize production and yield values.In 2019 and 2020, the temperature was above the regional average with occasional heatwaves; however, the precipitation exceeded 200 mm in the vegetation period of maize.For this study, three experimental fields were cultivated with maize on the Acqua Campus farm as shown in Figure 1.According to the CER data, the crop yield varied between 16.5 and 16.7 t/ha in 2018, 13.8 and 16.2 t/ha in 2019, and 18.8 and 19.9 t/ha in 2020, respectively.

Study Framework and Data
In this study, data from three maize fields in three years (2018,2019,2020) were used to set an advanced RS-based crop evapotranspiration estimation method to support a more accurate calculation of the water consumption of maize.The daily meteorological data (minimum and maximum air temperature, relative humidity, radiation, and wind speed) were used to estimate ET0 using the Penman-Monteith method.The crop coefficient was defined by developing linear and power regression equations between the NDVI, NDWI, NDRE LAI, and FAO Kc for the maize crop.Then VI-based Kc was used to calculate ETc, and the performances of the index-based crop evapotranspiration were compared to ETc derived by FAO Kc.
The analysis was conducted using a diverse set of datasets, incorporating Sentinel-2 satellite-derived NDVI, NDWI, NDRE, and LAI datasets, each possessing a high spatial resolution of 10 m.The LAI was calculated from the Sentinel-2 data using the methodology developed in the S2ToolBox within SNAP 5.0.S2ToolBox is an essential component of the SNAP 5.0 software and was specifically developed for LAI retrieval.The methodology uses the approach of an Artificial Neural Network (ANN) to derive these parameters from instantaneous observations of Sentinel-2.By utilizing a pre-trained neural network, the necessary biophysical variables are rapidly retrieved for each pixel within the selected Sentinel-2 image.The training database is generated using a Radiative Transfer Model (RTM) [57].
These datasets were acquired through the VULTUS API (https://api.vultus.se/graphql(accessed on 15. 02.2021)).The temporal resolution of the indices was set at five days, aligning with the satellite's revisit time.A comprehensive dataset, comprising 28 Vegetation Indices (VI) images and associated data, was annually collected, spanning from sowing to harvest.Only data free from cloud cover were considered for this study.This resulted in a final dataset of 996 observations.To ensure temporal coherence in the data, a Piecewise Cubic Hermite Interpolating Polynomial (PCHIP) approach was employed, resulting in a smoothed daily curve for indices throughout the vegetation period.While Acqua Campus of CER provided the weather data (daily minimum and maximum temperature, wind speed, relative humidity, and solar radiation) using an in situ weather station, as well as the Biologische Bundesanstalt, Bundessortenamt and CHemical industry (BBCH) [58] and LAI datasets (local observations).The weather data were used for calculating ET0, whilst BBCH and LAI were used to define the lengths of the phenological stages.Crop phenophases were assessed qualitatively by biweekly survey and sampling of maize plants in the field.The BBCH scale is used to identify the lengths of phenological development stages of plants [59].LAI measurement was carried out to assess the length of phenological stages through direct methods, which involves the removal of plants in

Study Framework and Data
In this study, data from three maize fields in three years (2018, 2019, 2020) were used to set an advanced RS-based crop evapotranspiration estimation method to support a more accurate calculation of the water consumption of maize.The daily meteorological data (minimum and maximum air temperature, relative humidity, radiation, and wind speed) were used to estimate ET 0 using the Penman-Monteith method.The crop coefficient was defined by developing linear and power regression equations between the NDVI, NDWI, NDRE LAI, and FAO K c for the maize crop.Then VI-based K c was used to calculate ET c , and the performances of the index-based crop evapotranspiration were compared to ET c derived by FAO K c .
The analysis was conducted using a diverse set of datasets, incorporating Sentinel-2 satellite-derived NDVI, NDWI, NDRE, and LAI datasets, each possessing a high spatial resolution of 10 m.The LAI was calculated from the Sentinel-2 data using the methodology developed in the S2ToolBox within SNAP 5.0.S2ToolBox is an essential component of the SNAP 5.0 software and was specifically developed for LAI retrieval.The methodology uses the approach of an Artificial Neural Network (ANN) to derive these parameters from instantaneous observations of Sentinel-2.By utilizing a pre-trained neural network, the necessary biophysical variables are rapidly retrieved for each pixel within the selected Sentinel-2 image.The training database is generated using a Radiative Transfer Model (RTM) [57].
These datasets were acquired through the VULTUS API (https://api.vultus.se/graphql(accessed on 15 February 2021)).The temporal resolution of the indices was set at five days, aligning with the satellite's revisit time.A comprehensive dataset, comprising 28 Vegetation Indices (VI) images and associated data, was annually collected, spanning from sowing to harvest.Only data free from cloud cover were considered for this study.This resulted in a final dataset of 996 observations.To ensure temporal coherence in the data, a Piecewise Cubic Hermite Interpolating Polynomial (PCHIP) approach was employed, resulting in a smoothed daily curve for indices throughout the vegetation period.While Acqua Campus of CER provided the weather data (daily minimum and maximum temperature, wind speed, relative humidity, and solar radiation) using an in situ weather station, as well as the Biologische Bundesanstalt, Bundessortenamt and CHemical industry (BBCH) [58] and LAI datasets (local observations).The weather data were used for calculating ET 0 , whilst BBCH and LAI were used to define the lengths of the phenological stages.Crop phenophases were assessed qualitatively by biweekly survey and sampling of maize plants in the field.The BBCH scale is used to identify the lengths of phenological development stages of plants [59].LAI measurement was carried out to assess the length of phenological stages through direct methods, which involves the removal of plants in the field on certain Plants 2024, 13, 1212 6 of 20 test areas and measuring their area in the laboratory.More details for LAI measurements are described by Nagy et al. [60].

VI-Based ET c Calculations
The crop evapotranspiration ET c was calculated using the following equation (ET c = ET 0 × K c ). ET 0 was defined by temperature, relative air humidity, solar radiation, wind speed, and air pressure according to the FAO Penman-Monteith equation [10] (Equation ( 1)). where: • ET 0 : reference crop evapotranspiration (mm day −1 ), Based on the BBCH and LAI, the lengths of the phenological stages were calculated for the three assessed years (Table 2).In general, the sowing time was in late March or early April, the harvesting was in September.The initial stage starts on the day of sowing and continues until about 10% of the soil cover is reached.The early root growth of maize makes it particularly sensitive to drought stress.In the early stages, evaporation is mainly in the form of soil evaporation and the leaf area is small.The mid-season is characterized by rapid growth of maize, from 10% soil cover to actual full soil cover.As the plant grows and covers a larger area of soil, transpiration gradually takes over as the primary mechanism, while evaporation becomes more limited, increasing water use.Water demand is highest in the mid-season stage due to high transpiration, increased biomass production, and grain formation, and K c reaches its maximum value.In the late season stage, transpiration-related water demand may decrease, but at this stage, drought stress can still have a significant negative impact on productivity [61][62][63].Since the determination method of K c presented in the FAO-56 document shows a trapezoid development, and the crop coefficient is defined as a constant throughout the initial and mid-stage, it does not reflect the natural continuous process of canopy development [32].This approach is referred to as "original" in this study.Therefore, in this study, the trapezoid character of the FAO K c was smoothed to have continuously varying characteristics by the application of two technical approaches curve with continuously varying values to be closer to natural crop development.One of the approaches was to cut 10 days around the breakpoints and the spline algorithm was implemented to estimate missing values.The solution is called "curved" in this study.The other approach was the so-called "MidPoint", where only the middle stage plateau was determined by its middle point and the length of the plateau phase.
To analyze the relationship between the FAO-56-based ("original", "MidPoint", "curved") K c values and the four Sentinel-2 derived indices (NDVI, NDRE, NDWI, and LAI), Pearson Product-Moment correlation was applied.Based on the results, FAO-56based K c was selected, which had the highest correlation with the Vis.Then linear and power regression models were developed based on VI data.Since certain Vis (such as NDVI) have proved to be an exponential nature compared to the other indices, both linear and power regression models were evaluated [61, [64][65][66].The use of standard linear or power regression models with standard estimation techniques is subject to several conditions regarding the explanatory (x) (which were the Vis) and output (y) variables (which were the K c values) and their relationship.
Using the Vis-based estimated K c and FAO-56 K c , standard ET c was calculated for all parcels for each year, resulting in biomass-derived evapotranspiration for the vegetation period of the maize fields.To evaluate the VI-based estimation ET c models, the following statistical indicators were selected: adjusted R 2 (R 2 adj ) (Equation ( 2)), Root Mean Square Error (RMSE) (Equation ( 3)), Normalized Root Mean Square Error (NRMSE) (Equation ( 4)), Mean Bias Errors (MBE) (Equation ( 5)) and Mean Absolute Error (MAE) (Equation ( 6)).RMSE is the sample standard deviation of the differences between predicted and actual values and NRMSE is useful to reflect the relative error between the modeled and measured soil moisture.The MBE was used to assess over or underestimation of VI-based ET c .The MAE evaluates the mean magnitude of the errors in predictions without considering their sign.The reason behind employing RMSE is its ability to measure accuracy to compare predicting errors of varying models for a particular dataset and not between datasets.The study used adjusted R 2 because it provides a measure of how well the predicted values are replicated by the model based on the variability from the actual values.All the parameters are defined as follows: ETc i A is the evapotranspiration based on theoretical K c values, ETc i P is the simulated evapotranspiration based on the spectral index derived and simulated K c values ETc − is the mean value, n is the total number of data points and k is the number of variables in the model.Adjusted R 2 : Root Mean Square Error: Normalized Root Mean Square Error: Plants 2024, 13, 1212 8 of 20 Mean Bias Error: Mean Absolute Error: In the initial phase (days 80-140), the mean NDVI values were low, with a mean value of 0.21.According to Lima et al. [70], NDVI values were low in the early part of the phenological cycle because of the higher radiation absorption in the near-infrared band during this period.In the crop development (days 141-179), NDVI values gradually increased and reached a maximum on day 174 (0.34).In the mid-season stage, values ranged from 0.56 to 0.67 until day 224 (late season).Reyes-González et al. [32] and Mebrie et al. [71] suggest that the crop coefficient varies constantly during the growing season and cannot that the leaves had little moisture content and that vegetation had not yet reached its full potential.In the crop development stage (days 141-179), the NDRE increased steadily as the vegetation grew, reaching a maximum on day 174 (0.82).During the mid-season (days 180-224), the mean NDRE values were similar to those around 0.76 at the crop development stage, similar to Liu et al. [69], who found, among other things, that the NDRE index has the highest correlation with canopy nitrogen concentration in summer maize.In the late-season stage (days 225-240), NDRE values showed a decreasing trend due to the gradual depletion of crop moisture.
In the initial phase (days 80-140), the mean NDVI values were low, with a mean value of 0.21.According to Lima et al. [70], NDVI values were low in the early part of the phenological cycle because of the higher radiation absorption in the near-infrared band during this period.In the crop development (days 141-179), NDVI values gradually increased and reached a maximum on day 174 (0.34).In the mid-season stage, values ranged from 0.56 to 0.67 until day 224 (late season).Reyes-González et al. [32] and Mebrie et al. [71] suggest that the crop coefficient varies constantly during the growing season and cannot remain constant throughout the mid-stage.In the late stage, the NDVI values showed a decreasing trend.Several studies have been conducted to investigate NDVI at different growth stages.Wang et al. [72] found similar values at different stages.The initial NDVI value of 0-0.2 increased with plant development up to 0.3-0.6,reaching a maximum of 0.7-0.8.Hatfield and Prueger [64] and Ji et al. [73] also divided the maize growing season into four growth phases and calculated phenological information metrics for each growth phase based on NDVI time series, while other authors have divided the growth phase into three, six, or even ten stages [74][75][76].
The trends and distribution of LAI values were similar to the previous indices.At the initial stage (days 80-140), the mean LAI values were low, around 0.5, as the vegetation was not fully developed, and the number of leaves was low.This statement aligns with the findings of Towers et al. [77], which supports the observation that at the initial stage, the mean LAI values were low due to underdeveloped vegetation and a low leaf count.In the crop development stage, the values increased as already described above and reached a maximum on day 174 (3.26).In the mid-season stage (days 180-224), the vegetation maintains a mean LAI value of around 2.7.The trend of data aligns with the findings of Zhang et al. [78] which supports the observation that LAI values follow an upward trend during the crop development stage and maintain a mean value during the mid-season stage in maize canopies.The LAI values also show a downward trend as leaf area decreases, like the other spectral indices.Similar trends in LAI values have been reported by several authors [79][80][81].
Overall, the trends and distributions of NDWI, NDRE, NDVI, and LAI values were similar, and the data followed well the vegetation development and environmental conditions.

Deriving VI-Based Crop Coefficient for the Vegetation Period of Maize
The lowest K c values were observed during the initial stage (0.4) based on FAO-56.In the case of the "curved" dataset, it typically starts to gradually increase from this initial value about 10 days earlier compared to the "original" and "MidPoint" values.However, the "original" K c values peak much earlier at the value of 1.2 in the crop development stage (days 166-170) than the other two datasets.Following the "original" K c values, the "curved" dataset reaches its maximum (1.2) approximately 10 days later, followed by the "MidPoint" dataset after 7-10 days.As illustrated in Figure 3, this plateau phase represents approximately a 50-day period in the original K c values, a 26-day period in the "curved" dataset for 2018 and 2019, and a 35-day period in 2020.In the case of the "MidPoint" dataset, this phase lasted for about only 8-9 days.stage (days 166-170) than the other two datasets.Following the "original" Kc values, the "curved" dataset reaches its maximum (1.2) approximately 10 days later, followed by the "MidPoint" dataset after 7-10 days.As illustrated in Figure 3, this plateau phase represents approximately a 50-day period in the original Kc values, a 26-day period in the "curved" dataset for 2018 and 2019, and a 35-day period in 2020.In the case of the "Mid-Point" dataset, this phase lasted for about only 8-9 days.There was a significant correlation between FAO-56-based Kc values and VIs, thus proper regression models were capable of describing the real crop evapotranspiration based on the evaluated vegetation index in this study (Table 3).The correlation between Vis and FAO-56-based Kc smoothed by the "MidPoint" method was found to be the highest, therefore "MidPoint" Kc was selected to use for setting regression models between the VIs and FAO-56 Kc.The correlation matrix indicates that LAI was less effective in reproducing the temporal pattern of the theoretical curves, while NDWI and NDVI seem to be the most effective in predicting Etc.Both linear and power regression models were set based on "MidPoint" Kc and the mean NDWI, NDRE, NDVI, and LAI values (Figure 4).The linear correlation between There was a significant correlation between FAO-56-based K c values and VIs, thus proper regression models were capable of describing the real crop evapotranspiration based on the evaluated vegetation index in this study (Table 3).The correlation between Vis and FAO-56-based K c smoothed by the "MidPoint" method was found to be the highest, therefore "MidPoint" K c was selected to use for setting regression models between the VIs and FAO-56 K c .The correlation matrix indicates that LAI was less effective in reproducing the temporal pattern of the theoretical curves, while NDWI and NDVI seem to be the most effective in predicting Etc.Both linear and power regression models were set based on "MidPoint" K c and the mean NDWI, NDRE, NDVI, and LAI values (Figure 4).The linear correlation between mean NDWI values and FAO crop coefficient of maize (R 2 ) was 0.83, indicating a strong positive correlation between the two variables.The R 2 value of the power trend line was even higher (0.88).However, we could only fit this to the positive integers, which cover the period from day 129 to day 234 of the year, while the total growth stage was from day 84 to day 249.The relationship between NDRE and FAO K c indicates a strong positive correlation for both linear and power fitting.The values were almost identical, R 2 being 0.8858 in the former and 0.8856 in the latter case.
The correlation between NDVI and FAO data for maize has been investigated in several studies.In this study, for the NDVI, almost the same can be stated as for the NDRE index, as the R 2 values were almost the same for both linear and power trend lines (R 2 = 0.886 and 0.8875) and show a strong correlation between NDVI and FAO K c .The relationship between the NDVI and the crop coefficient (K c ) for maize has been extensively studied.A strong linear association between the NDVI and FAO crop coefficient values was also discovered by Costa et al. [82]; the former's value was 0.86, while the latter was 0.79.Reyes-González et al. [66] discovered that the NDVI and K c FAO-56 had an even greater coefficient of determination of 0.97.Javed et al. [83] also reported that crop coefficients are highly correlated with satellite-derived NDVI values, suggesting a direct proportional relationship between the maize crop coefficient (FAO-56) and NDVI.Costa et al. [82] revealed a significant correlation between NDVI and the K c values reported in the FAO-56 report with R 2 = 0.794, an indication that there is a high correlation within these parameters.
Singh and Irmak [31] formulated a linear regression model to determine the correlation between NDVI and ET 0 -based K c .The results showed that the coefficient of variation for both NDVI and K c was lower for maize, soybeans, sorghum, and alfalfa at midseason compared to the early growing season as well as the late growth stage.Figure 5 illustrates the comparison between the FAO crop coefficient and crop coefficients estimated by linear and power functions for the growth stages.In the case of NDWI, in the initial stage (days 80-140), the linear calibrated Kc values were higher than the FAO Kc data, except on day 84, when the linear Kc was 0.379 and the FAO Kc was 0.4.The power-calibrated Kc data, as can be seen in Figure 5, could only be examined from day 129 onwards.In this short initial period (days 129-140), the power-calibrated Kc values and the FAO Kc values were almost the same.In the crop development stage (days 141-179), the power and linear calibrated Kc values were very similar (with a mean value of 0.99 in this stage) and had values higher than the FAO Kc (mean value of 0.94).In the mid-season stage (days 180-224), both FAO Kc and VI calibrated Kc values were above 1.Linear and power-calibrated values were again similar.The maximum crop coefficient value for the FAO data was on the 194th day (1.2), while for the power and linear calibrated data it was on the 174th day, 1.18 in the former case and 1.22 in the latter case.In the late season (days 225-250), values were again below 1, and similar to the mid-season, FAO Kc values were higher for both power and linear calibrated values.For the NDRE index, the linear and power-calibrated Kc values were almost equal in all stages.The FAO Kc data showed higher values than the calibrated data between days 91 and 115 of the initial stage (days 80-140), but for the remaining days of the initial stage and the whole crop development stage (days 141-179), the calibrated data were generally higher than the FAO Kc data.While in the initial stage, the mean of the calibrated data was 0.46 (for both power and linear), the FAO Kc was 0.41.The maximum values of power and linear calibrated Kc's were also the highest on the 174th day, 1.166 and 1.18, respectively.In the crop development stage, the mean of the calibrated Kc values for power and linear were The lowest R 2 was detected in the case of LAI between the index and the crop coefficient.The value for linear regression was 0.7712, whereas R 2 in power fitting was slightly lower at 0.7455.According to Kang et al. [84], the K c exhibited a rapid increase as the LAI approached 3 (linear R 2 = 0.9522), maintaining elevated values when LAI exceeded 2.5 (linear R 2 = 0.7612).The study by Jia and Wang [85] revealed a significant exponential relationship between LAI and crop coefficient K c (R 2 = 0.7055) under water deficit conditions, indicating the influence of vegetation indices on K c , and the impact of water stress on this relationship.
Figure 5 illustrates the comparison between the FAO crop coefficient and crop coefficients estimated by linear and power functions for the growth stages.In the case of NDWI, in the initial stage (days 80-140), the linear calibrated K c values were higher than the FAO K c data, except on day 84, when the linear K c was 0.379 and the FAO K c was 0.4.The power-calibrated K c data, as can be seen in Figure 5, could only be examined from day 129 onwards.In this short initial period (days 129-140), the power-calibrated K c values and the FAO K c values were almost the same.In the crop development stage (days 141-179), the power and linear calibrated K c values were very similar (with a mean value of 0.99 in this stage) and had values higher than the FAO K c (mean value of 0.94).In the mid-season stage (days 180-224), both FAO K c and VI calibrated K c values were above 1.Linear and power-calibrated values were again similar.The maximum crop coefficient value for the power and linear calibrated data were 1.19 in the former case and 1.29 in the latter case.If the mean values were taken into account, the mean of the FAO Kc values were higher (1.13 and 0.82) in the mid-season (180-224th day) and late-season stages (225-250th day).Among the calibrated crop coefficients, in the mid-season, the linear calibrated values were higher (with a mean value of 1.12, while the power was 1.08), while in the mid-season the power calibrated Kc values were slightly higher (0.65-with a mean value of 0.64).

Estimated Crop Evapotranspiration of Maize along the Growth Stages
The estimation of crop evapotranspiration for maize at different growth stages is a critical aspect of agricultural water management (Table 4).The minimum NDWI-based crop evaporation was on day 80 for both the power (1.34 mm/day ± 0.102) and linear (1.37 mm/day ± 0.049) calibrated data, while the maximum value for power was on day 185 (6.45 mm/day ± 0.367) and for the linear calibrated ET c data on day 170 (6.51 mm/day ± 0.434).
For NDRE, the minimum values were on day 95, which is still in the initial stage, with power at 1.36 mm/day ± 0.355 and linear calibrated ET c at 1.29 mm/day ± 0.293.The maximum ET c values were similar to the NDWI-based ET c , i.e., the power-calibrated ET c had its highest value at mid-season (6.76 mm/day ± 0.337 on day 185), while the linear calibrated ET c had its highest value towards the end of the crop development stage (6.4 ± 0.838 on day 170).
The evolution of NDVI-based K c values was very similar to NDRE-based K c values.The minimum value was on day 95, 1.31 mm/day ± 0.366 for the power calibration and 1.22 mm/day ± 0.286 for the linear calibration, and the maximum was on day 185 for the power calibrated K c values (6.53 mm/day ± 0.331) and on day 170 for the linear calibrated values (6.46 mm/day ± 0.669).
The power-calibrated LAI-based K c was lowest on day 80 (1.32 mm/day ± 0.015) and highest on day 185 (5.83 mm/day ± 0.492).In contrast, the linear calibrated K c values had a minimum on day 95 (1.27 mm/day ± 0.293) and a maximum on day 170 (5.93 mm/day ± 0.322).
In general, the vegetation index-based crop coefficient values were lowest in the initial stage, while the highest values were towards the end of the crop development stage (for linear calibrated, etc.) and in the mid-season stage (for power-calibrated ET c ).Based on the adjusted R-squared, the values closest to 1 were observed for the NDWI, NDRE, and NDVI indices, indicating that these models were able to explain a larger proportion of the variability in the data (Figure 6).For these three indices, the power regression model exceeded the linear regression model in general by a small margin.The LAI index had the lowest adjusted R-squared value, suggesting that the models could not fully account for the variability in the data.For LAI, the power regression model outperformed the linear regression model by a small margin; nonetheless, both models' predictions displayed very significant variability, indicating that the models are less accurate in predicting LAI compared to NDWI, NDRE, or NDVI.According to the RMSE, the NDWI, NDRE, and NDVI indices exhibit similar RMSE values across both regression models (ranging between 0.43 and 0.45 mm/day), indicating comparable performance for predicting these values (Figure 7).The LAI index demonstrated the highest RMSE values for both regression models.However, the linear regression model (0.5 mm/day) appears to provide a slight advantage compared to the power regression model (0.53 mm/day).
The results for NRMSE were similar to those for the statistical indices described so far (Figure 7).For the NDRE index, the NRMSE values were relatively close, with the linear regression model achieving slightly lower values compared to the power regression model (11.36%).This indicates that both models perform similarly in predicting NDRE values.The values for the NDVI index also showed a mixed pattern.The linear regression According to the RMSE, the NDWI, NDRE, and NDVI indices exhibit similar RMSE values across both regression models (ranging between 0.43 and 0.45 mm/day), indicating comparable performance for predicting these values (Figure 7).The LAI index demonstrated the highest RMSE values for both regression models.However, the linear regression model (0.5 mm/day) appears to provide a slight advantage compared to the power regression model (0.53 mm/day).For MAE, the lowest values were observed for NDVI (0.32 and 0.33), but NDWI and NDRE values were also close to this (0.34 and 0.35) (Figure 8).The standard deviation of the values ranged from 0.08 to 0.11.Overall, the NDWI, NDRE, and NDVI indices exhibited similar MAE values across both regression models, indicating comparable performance for predicting these values.The LAI index demonstrates the highest MAE values, for both regression models, with the power regression model achieving a mean value of 0.40 mm/day and the linear regression model achieving a mean value of 0.37 mm/day.
In the case of MBE (Figure 8), the estimates of the linear regression model show mean smaller deviations from the actual values than the estimates of the power regression model.This was particularly the case for NDWI and LAI, where the mean NDWI for the power regression model was −0.04 mm and −0.005 mm for the linear model, while for LAI it was −0.042 mm and −0.003 mm.The NDRE and NDVI indices were closer to 0 for both the power and linear regression models.The results for NRMSE were similar to those for the statistical indices described so far (Figure 7).For the NDRE index, the NRMSE values were relatively close, with the linear regression model achieving slightly lower values compared to the power regression model (11.36%).This indicates that both models perform similarly in predicting NDRE values.The values for the NDVI index also showed a mixed pattern.The linear regression model achieves relatively lower mean NRMSE values (11.56%) compared to the power regression model (11.48%), suggesting a small advantage for the linear regression model in predicting NDVI values.The LAI index exhibited the highest NRMSE values for both regression models.
For MAE, the lowest values were observed for NDVI (0.32 and 0.33), but NDWI and NDRE values were also close to this (0.34 and 0.35) (Figure 8).The standard deviation of the values ranged from 0.08 to 0.11.Overall, the NDWI, NDRE, and NDVI indices exhibited similar MAE values across both regression models, indicating comparable performance for predicting these values.The LAI index demonstrates the highest MAE values, for both regression models, with the power regression model achieving a mean value of 0.40 mm/day and the linear regression model achieving a mean value of 0.37 mm/day.
In the case of MBE (Figure 8), the estimates of the linear regression model show mean smaller deviations from the actual values than the estimates of the power regression model.This was particularly the case for NDWI and LAI, where the mean NDWI for the power regression model was −0.04 mm and −0.005 mm for the linear model, while for LAI it was −0.042 mm and −0.003 mm.The NDRE and NDVI indices were closer to 0 for both the power and linear regression models.

Conclusions
In the Mediterranean, water competition among different stakeholders becomes more intense and affects agricultural production, especially in areas of high-intensive use like northern Italy.Changes in temperature patterns, rainfall, and extreme weather events affect agricultural outputs and make it difficult to plan optimal distribution of water and could lead to food insecurity.Effective water management becomes critical to addressing these challenges.For solving these issues, modern technological approaches such as RS provide effective solutions.In this study, researchers have created an innovative RS-based approach for calculating crop coefficients and evapotranspiration based on a case study site in Emilia Romagna, Italy.Through the use of remotely sensed indices coupled with climatic data, the accuracy in estimating crop evapotranspiration was substantially increased for maize fields.
The comparative study of NDVI, NDWI, NDRE, and LAI values in different stages of maize growth helps understand the changes that occur during vegetation development under changing environmental conditions.In addition, the relationship between these indices and Kc reveals their suitability in predicting crop evapotranspiration.The observed strong positive correlations between NDVI, NDWI, NDRE, and LAI values with FAO Kc indicate the possibility of using RS data to obtain reasonable estimates for crop water requirements.Regression models of NDVI, NDWI, and NDRE exhibit slightly different predictive performances in determining Kc values with power regression generally being superior to linear regressions.The most significant correlations are found for NDVI and NDWI, indicating their strength in predicting Kc values in maize.Nevertheless, there are

Conclusions
In the Mediterranean, water competition among different stakeholders becomes more intense and affects agricultural production, especially in areas of high-intensive use like northern Italy.Changes in temperature patterns, rainfall, and extreme weather events affect agricultural outputs and make it difficult to plan optimal distribution of water and could lead to food insecurity.Effective water management becomes critical to addressing these challenges.For solving these issues, modern technological approaches such as RS provide effective solutions.In this study, researchers have created an innovative RS-based approach for calculating crop coefficients and evapotranspiration based on a case study site in Emilia Romagna, Italy.Through the use of remotely sensed indices coupled with climatic data, the accuracy in estimating crop evapotranspiration was substantially increased for maize fields.
The comparative study of NDVI, NDWI, NDRE, and LAI values in different stages of maize growth helps understand the changes that occur during vegetation development under changing environmental conditions.In addition, the relationship between these indices and K c reveals their suitability in predicting crop evapotranspiration.The observed strong positive correlations between NDVI, NDWI, NDRE, and LAI values with FAO K c indicate the possibility of using RS data to obtain reasonable estimates for crop water requirements.Regression models of NDVI, NDWI, and NDRE exhibit slightly different predictive performances in determining K c values with power regression generally being superior to linear regressions.The most significant correlations are found for NDVI and NDWI, indicating their strength in predicting K c values in maize.Nevertheless, there are lower correlation values for LAI which indicates that the predictive accuracy of this index is less than other indices.
Generally, NDWI, NDRE, and NDVI-based ETc estimations outperform consistently compared to LAI with power regression providing better results.In the end, it can be concluded that VI (NDWI, NDRE, NDVI, and LAI)-based crop evapotranspiration ETc are indicating different performances in predicting water requirements for crops compared to FAO-56 Kc-based ETc.NDWI, NDRE, and NDVI have good explanatory power with relatively accurate predictions for both models of regression while the LAI index shows some limitations in its accuracy.The power regression model usually slightly exceeds the linear regression model, NDWI, NDRE, and NDVI have better predictive accuracy than LAI.These results emphasize the importance of vegetation indices, especially NDWI and NDVI as promising approaches to estimate crop water requirements with further improvements required for LAI-based predictions.
In conclusion, the analysis shows that NDWI and NDVI indices are effective in monitoring maize growth stages and estimating crop evapotranspiration.Using RS data and regression modeling techniques provides appropriate tools for the improvement of agricultural water management practices to achieve higher crop productivity.More research and improvement of these methodologies can lead to more accurate and consistent strategies for crop water management in agricultural systems.

Figure 1 .
Figure 1.Experimental fields on the Acqua Campus farm.

Figure 1 .
Figure 1.Experimental fields on the Acqua Campus farm.

Figure 2 .
Figure 2. The mean NDWI, NDRE, NDVI, and LAI values over the three years (2018-2020).The trends and distribution of NDRE values are similar to NDWI values.The mean NDRE values during the initial phase (days 80-140) were relatively low, at 0.2.This suggests

Figure 5 .
Figure 5.Comparison of FAO-56 crop coefficient, and linear and power-calibrated crop coefficient of maize with NDWI, NDRE, NDVI, and LAI indices.

Figure 5 .
Figure 5.Comparison of FAO-56 crop coefficient, and linear and power-calibrated crop coefficient of maize with NDWI, NDRE, NDVI, and LAI indices.

Plants 2024 ,
13,  x FOR PEER REVIEW 15 of 21 regression model exceeded the linear regression model in general by a small margin.The LAI index had the lowest adjusted R-squared value, suggesting that the models could not fully account for the variability in the data.For LAI, the power regression model outperformed the linear regression model by a small margin; nonetheless, both models' predictions displayed very significant variability, indicating that the models are less accurate in predicting LAI compared to NDWI, NDRE, or NDVI.

Figure 6 .
Figure 6.Adjusted R 2 between ETc estimated with the vegetation indices (NDWI, NDRE, NDVI, and LAI) and that estimated with FAO-56 K c .

Figure 8 .
Figure 8. MAE and MBE between ETc estimated with the vegetation indices (NDWI, NDRE, NDVI, and LAI) and that estimated with FAO-56 Kc.

Figure 8 .
Figure 8. MAE and MBE between ETc estimated with the vegetation indices (NDWI, NDRE, NDVI, and LAI) and that estimated with FAO-56 Kc.

Table 1 .
Summary of mean values for the soil parameters.

Table 2 .
The lengths of the phenological stages.

Table 3 .
Correlation coefficients between estimated Kc curves.

Table 3 .
Correlation coefficients between estimated K c curves.